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Arrangement 



for suppressing an interfering component of an input signal* 



10 



The present invention is related to an arrangement for suppressing an 
interfering component in an input signal, said arrangement comprises means for deriving an 
estimate of the interfering component, subtraction means for determining a difference s.gna. 
from the input signal and the estimate of the interfering component. 

The present invention is also related to a loudspeaking telephone . 
An arrangement according to the preamble is known from US 5,390,250. 
This arrangement has widespread applications such as acoustic echo 
cancellers and noise cancellers. In these application there is an input signal in which an 

interfering component is present. 

This interfering component can e.g. be a signal received from the far end 
in a loud speaking telephone system. This far end signal is reproduced by a loudspeaker and 
is received by a microphone together with the near end signal e.g. from a local speaker. The 
input signal is now a signal derived from the output signal of the microphone. Because 
amplifiers are included in the transmit path, the «oo P gain for a specific frequency may be 
15 greater than 1, resulting in oscillations. If the acoustic feedback is smaller than 1. oscillation 
will not occur, but after a certain delay an echo of the signal applied to the input of the send 
path will appear at the output of the receive path via the far-end echo path. In telephony th.s 
means that a speaker hears his own voice delayed by a specific period of tune. Tins 
phenomenon is experienced as extremely annoying especially in case of long delays. 
20 to prevent this undesired feedback, in an acoustic echo canceller a rephca 

of the undesired component is derived from the far end signal by means of an adaptive filter. 
Said replica is subtracted from the input signal, in order to eliminate the undes.red 
component signal. 

In case of a noise cancelling system the undesired component is a noise 
25 signal coming from a noise source, e.g. the noise of a running engine in a car. To eliminate 
the noise signal, a replica of the noise signal is derived from a reference noise signal 
received from a reference transducer by means of an adaptive filter. Again, this replica ,s 

subtracted from the input signal. 

A problem with adaptive filters is the limited ability to track fast changes 
30 of the transfer function to be reproduced. Such a fast change can occur due to the movement 
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of a person in the room in which the echo canceller or noise canceller is used. This can 
result into a reduced amount of suppression of the undesired component, a problem that ,s 

not satisfactorily solved up to now. 

The object of the present invention is to provide an arrangement accord.ng 

5 to the preamble in which the ability to cope with changes in the transfer function to be 
reproduced by the adaptive filter has been drastically been improved. 

Therefor the arrangement according to the invention is characterised in 
that the arrangement comprises means for determining an estimate for the frequency 
spectrum of the interfering component, and in that the arrangement comprises a filter having 

10 a reduced transfer function for at least one frequency range in dependence on the frequency 

spectrum of the interfering component. 

By determining the frequency spectrum of the interfering component and 
by the use of a filter for selectively attenuating the frequency ranges corresponding to the 
frequency determined, an additional suppression of the interfering component is obtained. 
15 Experiments have shown surprisingly that the use of the above mentioned filter has almost no 
perceptual effect on the desired component of the input signal. Also these experiments 
revealed that the additional suppression of the undesired component, allows very rapid 
changes in the transfer function to be reproduced without hardly any perceptual effect. 

An embodiment of the invention is characterised in that the filter is 
20 arranged for deriving an output signal from the difference signal. 

If the filter is arranged for deriving the output signal in dependence on the 
difference signal, the operation of the means for deriving the replica of the undesired 
component is fully decoupled from the filtering operation. This has as advantage that the 
convergence properties of the means for deriving the replica of the undesired component are 

25 independent of the additional filtering. 

A further embodiment of the invention is characterised in that the means 
for determining an estimate for the frequency spectrum of the interfering component are 
arranged for deriving the estimate for the frequency spectrum of the interfering component 
from the estimate of the interfering component. 
3Q At the output of the means for estimating the interfering component 

generally a good estimate of said interfering component is present. Using this signal for the 
determination of the frequency spectrum, results in a reliable estimation of said frequency 

SpeCtrUm ' A fufther embodim ent of the invention is characterised in that the filter 
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has a transfer function having an absolute value dependent on the difference of the amplitude 
of the frequency of the frequency spectrum of the input signal and the amplitude of the 
frequency spectrum of the estimate of the interfering component. 

An easy way of deriving the transfer function of the filter is the 
5 subtraction of the estimate of the amplitude spectrum of the undesired component from the 
amplitude spectrum of the input signal. This subtraction operation is normally performed in 
the frequency domain, but the present invention is not limited thereto. 

A very simple embodiment of the invention is characterised in that the 
means for determining an estimate for the frequency spectrum of the interfering component 
10 are arranged for determining the frequency of at least one peak in the frequency spectrum, 
and in that the filter is arranged for attenuating components in a frequency range around the 

frequency of the at least one peak. 

The frequency peak or peaks can e.g. be determined by calculation of 
linear prediction parameters. The filter can have a transfer function being dependent on the 

15 prediction parameters. 

The invention will now be explained with reference to the drawings. 

Herein shows: 

Fig. 1. a first embodiment of an arrangement according to the invention 

used in an echo canceller; 
20 pig. 2, a second embodiment of an arrangement according to the 

invention used in a noise canceller; 

Fig. 3, the position of the microphones to be used in the embodiment 

according to Fig. 2; 

Fig. 4, an embodiment of the means for deriving a replica of the 

25 interfering signal to be used with the embodiment according to Fig. 1 or Fig. 2; 

Fig. 5, a first embodiment of the combination of the means for estimating 
the frequency spectrum of the interfering component and the filter; 

Fig. 6, a second embodiment of the combination of the means for 
estimating the frequency spectrum of the interfering component and the filter. 

30 In the echo canceller according to Fig. 1. a signal xfk] is connected to a 

first input of a signal processing arrangement 6. In the arrangement 6. the signal x[k] is 
applied to an input of an amplifier 2 and to an input of the means for deriving an estimate of 
the interfering component, said means being here an adaptive filter 10. The output of the 
amplifier 2 is connected to an input of a loudspeaker 4. 
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The output of a microphone, carrying the input signal z[k] of the arrange- 
ment is coupled to a first input of a subtracter 14. The output of the adaptive filter 10 is 
connected to a second input of the subtracter 14, and to an input of the means for 
determining an estimate of the frequency spectrum of the interfering signal, further to be 
5 referred to as spectrum determining means 12. The output of the subtracter 14 is coupled to 
an input of a filter 16. The output of the spectrum determination means 12 is connecter to a 
control input of the filter 16. At the output of the filter 16 an estimate r'[k] of the signal s|k1 
is available. 

In the echo canceller according to Fig. 1, the interfering component in the 
10 input signal z[k] is the signal e[k]. This signal e[k] is an echo signal caused by the 

reproduction of the far end signal by the loudspeaker 4. The adaptive filter 10 is arranged for 
deriving a replica e[k] of the signal e[k]. This is in general done by choosing the coefficients 
of the adaptive filter for minimising the correlation between the signal rfk] and the signal 
x[k] There exist several well known time domain algorithms for adjusting the coefficients of 
15 the adaptive filter, such as the LMS (Least Mean Square) algorithm, the NLMS (Normalised 
Least Mean Square) algorithm and the RLS (Recursive Least Square) algorithm. It is also 
possible that the adaptive filter operates in the frequency domain. 

The spectrum determining means 12 determine the frequency spectrum of 
the output signal of the adaptive filter 10. From this frequency spectrum, the setting for the 
filter 16 is determined. The filter 16 is arranged for suppressing the spectral components of 
the output signal of the subtracter having a strong contribution from the interfering s.gnal. It 
is conceivable that the input of the spectrum determining means is connected to the mput of 
the adaptive filter 10 instead to its output, because the signal x[k) is also a reasonable 
estimate for the interfering component in the input signal. It is also possible that the estimate 
25 for the interfering component is made during the absence of the signal s[k]. In this case a 
speech activity detector has to be added to the arrangement. With respect to the position of 
the filter 16 it is observed that it is possible that the filter 16 is present between the m.cro- 

phone 8 and the subtracter 14. 

In the noise canceller according to Fig. 2, the output of a microphone 20, 
30 carrying output signal z W , is connected to a first input of a signa. processing arrangement 
28. The signal z[k] is the input signal of the arrangement. In the arrangement 28 the s.gnal 
z[k] is applied to a' first input of a subtracter 38. 

An output of a microphone 22. carrying output signal x.M is connected 
to an input of an adaptive filter 30. The output of the adaptive filter 30 is connected to a first 
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input of an adder 36. An output of a microphone 24, carrying output signal x 2 [k] is 
connected to an input of an adaptive filter 32. The output of the adaptive filter 32 is 
connected to a second input of the adder 36. An output of a microphone 26, carrying output 
signal x 3 [k] is connected to an input of an adaptive filter 34. The output of the adaptive filter 
5 34 is connected to a third input of the adder 36. 

The output of the adder 36 is connected to a second input of the subtracter 
38 and to an input of spectrum determining means 42. The output of the subtracter 38. 
carrying output signal r[k] is connected to an input of a filter 40. an output of the spectrum 
determination means 24 is connected to a control input of the filter 40. At the output of the 
10 filter the signal r'[k] substantially free from the interfering component is available. 

The noise canceller according to Fig. 2 comprises four microphones 20. 
22, 24 and 26. The microphone 20 is the primary microphone intended for receiving signals 
from a desired speaker. The microphones 22, 24 and 26 may be assigned to competing 
speakers whose speech is not to be passed to the output, but the microphones 22, 24 and 26 
15 may also be reference microphones for receiving the interfering signal at three different 
positions in space. It goes without saying that the present invention is not limited to the 
application of three microphones for receiving the interfering signal. One microphone for 
receiving the interfering signal can be sufficient, but adding more microphones can improve 
the performance of the noise canceller. Each of the microphones 22, 24 and 26 is followed 
20 by an adaptive filter 30, 32 and 34 respectively.The combined signal ylk] derived from the 
output signals from the three adaptive filters 30, 32 and 34 is subtracted from the signal z[kj 
by the subtracter 38. The adaptive filters 30, 32 and 34 are individually adapted in order to 
minimise the correlation between the output signal of the subtracter 38 and its respective 
input signal. The same algorithms as in the adaptive filter 10 in Fig. 1 can be used. The 
25 spectrum determining means 42 determine the frequency spectrum of the output signal of the 
adder 36. The determined frequency spectrum is used to control the filter 40 for additional 
suppression of the interfering component in the output signal of the subtracter 38. 

Fig. 3 shows a possible arrangement of the four microphones 20, 22, 24 
and 26 for use in an arrangement according to Fig. 2. This arrangement is suitable for a 
30 plurality of speakers. Due to the directivity of each of the microphones, the microphones 22, 
24, 26 only receive reflections of the signal s instead of the signal s itself. This already leads 
to a suppression of the interfering component in the output signal of the noise canceller. By 
using the adaptive filters 30, 32 and 34 and the filter 40, the interfering components from the 
other speakers is effectively suppressed as is explained above. 
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In the implementation of the adaptive filter according to Fig. 4, the signal 
xjtk] is applied to an input of a frequency domain adaptive filter 50, and to an input of a 
time domain programmable filter. The input signal z[k] is applied to a first input of a 
subtracter 58 and to a first input of a subtracter 56. The output of the time domain 
5 programmable filter 54 is connected to a second input of the subtracter 58. At the output of 
the subtracter, the output signal is available. 

An output of the frequency domain adaptive filter 50 is connected to a 
second input of the subtracter 56. The output of the subtracter 56 is connected to a second 
input of the frequency domain adaptive filter 50. A second output of the frequency domain 
10 adaptive filter is connected to an input of an Inverse Fast Fourier Transformer 52. Said 
output of the frequency domain adaptive filter 50 carries the parameters of said frequency 
domain adaptive filter 50. The output of the Inverse Fast Fourier Transformer is connected 
to an input of the time domain programmable filter 54. 

The adaptive filter 50 is realised in the frequency domain because a time- 
15 domain adaptive filter for the necessary length of the impulse response is of much greater 
complexity than a frequency-domain adaptive filter. In addition, a time-domain adaptive filter 
has the disadvantage that adapting its impulse response of the echo path for signals having a 
strong auto correlation takes much longer than for a frequency-domain adaptive filter for 
which a decorrelation that is simple to implement is used, introducing such decorrelation into 
20 a time-domain adaptive filter would involve much greater complexity. 

In the frequency domain filter 50, blocks of samples of the signals Xj[k] 
and r"[kl are subjected to an orthogonal transform e.g. a Fast Fourier Transform to obtain 
the signals Xjlp.m] and R"[p;m], in which p is the normalised frequency, and m is the block 
number. In the frequency-domain filter 50 the correlation measure X^p^m] - R"[p;m] is 
25 calculated. Subsequently the frequency-domain coefficients Wlp;mJ of the filter 50 are 
updated according to the LMS algorithm. For the frequency-domain coefficients can be 
written: 



W[ P; m]=Wtp;nJ-l]+MtP;«'] ' x i [ P !in] R " [p;tn} 



(1) 



30 



In (1) fx is a gain factor which may be dependent on the power of Xilp;m]. The actual 
filtering takes place by multiplying the signal X^m] by the coefficients WLp;m] 
Subsequently the filtered signal is subjected to an Inverse Fast fourier transform to obtain the 
time-domain filtered signal to be applied to the second input of the subtracter 56. The 
frequency domain coefficients from the adaptive filter 50 are transformed into time-domain 



WO 97/45995 7 PCT/IB97/00564 

filter coefficients. These time domain filter coefficients are applied to the time-domain 
programmable filter 54, which actually determines the replica of the interfering component. 
The filter 9 is realised in the time domain because the signal to be filtered is not subjected to 
an additional delay in contradistinction to a frequency domain adaptive filter in which a 
5 certain additional delay is unavoidable. 

In the combination of the means 12 for estimating the frequency spectrum 
of the interfering component and the filter 16 according to Fig. 5, the signal y[k] is applied 
to a window processor 60. The output of the window processor 60 is connected to an input 
of an Fast Fourier Transformer 62. The combination of the window processor 60 and the 
10 Fast Fourier transformer 62 constitutes the means for determining an estimate of the 

frequency spectrum of the interfering component.The output of the Fast Fourier Transformer 
is connected to an amplitude detector 64. the output of the amplitude detector 64 is connected 
to a first input of a calculator 78. 

The signal y[k] is applied to an input of a window processor 66. The 
15 output of the window processor 66 is connected to an input of a Fast Fourier transformer 68. 
The output of the Fast Fourier Transformer 68 is connected to an input of an amplitude 
detector 74 and to a second input of the calculator 78. The output of the amplitude detector 
74 is connected to a third output of the calculator 78. 

The signal z[k] is applied to a window processor 70. The output of the 
20 window processor 70 is connected to an input of an Fast Fourier Transformer 72. The output 
of the Fast Fourier Transformer 72 is connected to an amplitude detector 76. The output of 
the amplitude detector 76 is connected to a fourth input of the calculator 78. 

The output of the calculator 78 is connected to an input of an Inverse Fast 
Fourier transformer 80. At the output of the Inverse Fourier Transformer 80 the output 

25 signal r*[k] is available. 

In the window processor 66, the most recent block of 128 samples of the 
signal r[k] is placed together with the most recent block but one of 128 samples into a buffer 
memory of 256 samples. The content of the buffer memory is windowed with a so-called 
Hanning window. Subsequently the most recent block of 128 samples of the signal r|k] is 

30 stored for used in the next window operation, and the most recent block but one of 128 
samples is discarded. The operation of the window processors 60 and 70 is the same as the 
operation of the window processor 66. 

The Fast Fourier Transformers 62, 68 and 72 calculate a 256 points FFT 
from their respective input signals to obtain the signals Y(«), R(w) and Z(u). The amplitude 
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detectors 64, 74 and 76 determine the amplitude of the frequency spectrum according to: 

R a iu,) = |»{K<«>} I + |9f{R(«))l (3) 

Z a (to) = |»{Z(«)}| + |9{Z(«))I (4) 

In (2), (3) and (4) W{x} indicates the real part of x, and 3{x} indicates the imaginary part of 
x. The approximation of the amplitude of the frequency spectra according to (2). (3) and (4) 

is used to simplify the calculations. 

The calculator 78 starts with calculating a first estimate of the amplitude 
of the frequency spectrum of the signal S. A possible expression for the amplitude spectrum 
S a of S is: 



[o, (z a -yW»))<f\ 



< 5 > 

S fl (w) =MAX 



10 



In (5) a and B are constants normally having a value between 0.5 and 2. y is a constant 
having a value around 1 or larger than 1. A suitable value is 1.5. The MAX function is 
introduced in (5) for preventing the amplitude spectrum to assume negative values. Expen- 
ments have shown that simplification of (5) by choosing a and B equal to 1 does not result ,n 
any audible degradation. Consequently (5) changes into: 



S„(w) 



,„ order to reduce the filtering action according to (6) in cases the adaptive filter 10 provides 
sufficient suppression of the interfering component, a spectra, floor is introduced in (6) 
15 according to: 

Vw: S a (w) = MAX[S a {0!) ,C Z a (u) } (7) 

In (7) C is a constant having a value around 0.003. 

Because the amplitude spectrum VM-SM of the filtered signal r'M is equal to 
H( W ).R». for Hfc.) can be written: HM- W^). ■« — » ,imit 
value of H(u), HM is made equal to MIN(1,S»/R»). Before H(«) is used to filter the 
20 signal r[k], a smoothed version of H(») is derived according to: 
In (8) m is the block number and » is a constant with value of 0.7. 



«> I 
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The frequency spectrum of the signal r'[k] is now determined from: 

r' (w) = R(cu) H m (u>) (9) 

The output signal of the calculator 78 is transformed into the time domain by the 256 points 
Inverse Fourier transform 80. The first 128 points of the Inverse FFT are added to the final 
128 points of the Inverse FFT previously calculated. The final 128 points of the present 
5 Inverse FFT are stored for combination with the first 128 points of the next IFFT to be 
calculated. 

In the arrangement 5. the signal Y» is determined directly from the 
output signal y[k] of the adaptive filter 10. However it is possible to obtain an even better 
estimation of Y». Normally the length of the impulse response of the adaptive filter will 
10 be limited. Consequently the echo signal present in the signal z[k) will show a tail being not 
reproduced by the adaptive filter. However, it is possible to include this tail in the estimation 
of Y a (w) as will be explained below. 

For the amplitude spectrum Y» of the echo signal can be written: 

Y a l") - Y af (w) * Y ">" (w) (10> 

In (10), Y af (u>) is the amplitude spectrum of the output signal of the adaptive filter, and 
15 Y^M is die amplitude spectrum of the residual echo. If it is assumed that the impulse 
response of the echo path decreases exponentially, the amplitude spectrum of the residual 
echo signal can be estimated from the output signal of the adaptive filter according to: 

Ytail [»]<«> =<* [»-lH«*> ^'V^P 1 '"' (11> 

In (11) cc and B are constants, m is an index indicating subsequently determined amplitude 
spectra Y af (a>), and p is the length of the adaptive filter expressed in a number of update 
20 periods of Y>>). In case the signal Y is stationary, (11) can be approximated by 

tm] (co) = a -Yuui [Jn-11 <«> + & *<tf™ <<•>> (12) 



*tail 



For a filter having an impulse response of 1024 samples, and if an update period of 128 
samples is used for Y>», suitable values for a and 11 are 0.79 and 0.04 respectively. 

In the combination of the means 12 for estimating the frequency spectrum 
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of the interfering component and the filter 16, the signal y[k] is applied to a linear predictive 
analyser 84. This linear predictive analyser determines a plurality of prediction coefficients 
art]. The transformer arrangement 86 performs a transformation to the prediction coefficient 
such as a bandwidth expansion according to a'[i]=a[i] • 8 l in which 6 is a number smaller 
than 1. The transformed prediction coefficients are used to control the filter 16, in order to 
suppress the frequency components for which the interfering component has a substantial 
amount of energy. A suitable filter could have the transfer function: 



l-£p f ati] z" 1 ' 
I °i 

- N 



»!«)- ^ 5 (13) 

z- ) 



i =1 



In (10) p. q and ft are constants smaller than 1 and p being smaller than q. (10) is based on 
the inverse transfer function of a post filter used in speech coding. The object of such a post 
filter is to enhance the part of the spectrum in which speech components are present, and to 
suppress noise components in the part of the spectrum in which no speech components are 
present. By taking a filter with an inverse transfer function, the part of the spectrum in which 
no interfering component is present is enhanced, and the part of the spectrum in which a 
strong interfering component is present is attenuated. 
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CLAIMS 



1 Arrangement for suppressing an interfering component in an input signal, 
said arrangement comprises means for deriving an estimate of the interfering component, 
subtraction means for determining a difference signal from the input signal and the estimate 
of the interfering component, characterised in that the arrangement comprises means for 
determining an estimate for the frequency spectrum of the interfering component, and in that 
the arrangement comprises a filter having a reduced transfer function for at least one 
frequency range in dependence on the frequency spectrum of the interfering component. 

2 Arrangement according to claim 1. characterised in that the filter is 
arranged for deriving an output signal from the difference signal. 

3 ' Arrangement according to claim 1 or 2, characterised in that the means 
for determining an estimate for the frequency spectrum of the interfering component are 
arranged for deriving the estimate for the frequency spectrum of the interfering component 
from the estimate of the interfering component. 

4 Arrangement according to one of the previous claims, characterised in that 
the filter has a transfer function having an absolute value dependent on the difference of the 
amplitude of the frequency of the frequency spectrum of the input signal and the amplitude of 
the frequency spectrum of the estimate of the interfering component. 

5 . Arrangement according to one of the claims 1 ,2 or 3, characterised in that 

the means for determining an estimate for the frequency spectrum of the interfering 
) component are arranged for determining the frequency of at least one peak in the frequency 
spectrum, and in that the filter is arranged for attenuating components in a frequency range 
around the frequency of the at least one peak. 

6 Loudspeaking telephone comprising an arrangement according to one of 

the claims 1,2,3,4 or 5. 
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